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Large prebiotic molecules in space: photo-physics of acetic acid and 
its isomers 
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ABSTRACT 

An increasing number of large molecules have been positively identified in space. Many 
of these molecules are of biological interest and thus provide insight into prebiotic organic 
chemistry in the protoplanetary nebula. Among these molecules, acetic acid is of particular 
importance due to its structural proximity to glycine, the simplest amino acid. We compute 
electronic and vibrational properties of acetic acid and its isomers, methyl formate and glyco- 
laldehyde, using density functional theory. From computed photo-absorption cross-sections, 
we obtain the corresponding photo-absorption rates for solar radiation at 1 AU and find them 
in good agreement with previous estimates. We also discuss glycolaldehyde diffuse emission 
in Sgr B2(N), as opposite to emissions from methyl formate and acetic acid that appear to be 
concentrate in the compact region Sgr B2(N-LMH). 
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1 INTRODUCTION 

The richness of interstellar chemistry has been growing steadily 
with the variety of objects and regions observed. The excita- 
tion and abundances of molecules contain key information on 
the physical structure and evolution of the host regions. Through 
molecules, we can trace the cycle of matter for interstellar space 
into stars and planets, and back again into the interstellar space 
dHartquist & Williams! 199^ . 

There is an increasing body of evidence for the existence of 
large molecules in the interstellar medium and in the interplanetary 
space. Interferometric observations of high-mass star forming re- 
gions in molecular clouds have revealed hot molecular cores, short- 
lived remnants of clouds not incorporated into the newly born mas- 
sive stars. These hot cores contain within them the evaporated ma- 
terial of ices deposited on dust grain surfaces during the collapse, 
and observations show a very interesting chemistry. In particular, 
there exist substantial column densities of large partly hydrogen- 
saturated molecules, many of them being of pre-biotic interest. 
Some of these species have been also obser ved in comets and em - 
bedded in minor bodies of the solar system dChvb a & Hand 20051) . 
The study of biologically interesting large interstellar molecules of- 
fers the exciting opportunity of learning more about the chemical 
evolution preceding the onset of life on the early Earth 4.5 billion 
years ago. Comets may be important carriers of prebiotic chem- 
istry, and relevant agents in the delivery of complex organics to 
early Earth, as well as to newly formed planets. 
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The Sgr B2 molecular cloud complex is the prime target in 
the search for complex species, in particular in a hot core, the 
so-called Large Molecule Heimat, Sgr B2(N-LMH), within the 
more extended molecular cloud. In this compact source, smaller 
than the P ort cloud (~ 0.08 pc) with a mass of several thou- 
sands M© dMiao&Snvdeilll997t) , an extraordinary number of 
complex organics have been observe d to exhibit very hi gh col- 
umn densities (e.g. amino acetonitrile, lBellocheetafl2 008). Large 
partly hydrogen-saturated species challenge the completeness of 
the standard ion-neutral scheme in interstellar chemistry, suggest- 
ing that reactions on du st grains are involved in their formation (e.g. 
iBennett & Kaisej|2007h . 

Of the chemical species detected so far, particular attention 
has been paid to the formation of different isomer groups. In 
this work, we focus on C2H4O2, i.e. acetic acid (CH3COOH), 
glycolaldehyde (HCOCH 2 OH), and methyl formate (HCOOCH3), 
because of their potential role in the origin of life (e.g., 
IWachtershause^l200d : lchvba & Handll2005h . Glycolaldehyde, the 
simplest of the m onosaccharide sugars, has first been detected by 
iHollis et all d2000h towards Sgr B2(N-LMH), and its most recently 
determined column density in that source is 5.9 x 10 13 cm~ 2 
jHalfen et alj 12006 ). It has bee n rece ntly observed outside the 
galactic center by Beltran et al.l j2009l) towards the hot molecu- 
lar core G31. 41+0.31, with the emission coming fr om the hot and 
dense regio n closest to the protos tars. In addition, ICrovisier et al.l 
d2004l) and iDespois et al.l d2005l) presented an upper limit for 
glycolaldehyde abundance in the comet Hale-Bopp. Acetic acid 
shares the C— C— O backbone with glycine, fro m which it dif- 
fers by an amino group (NH2). First detected by iMehringer et ail 
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j 19971) in Sgr B2(N-LMH), acetic acid shows a colum n density 
of 6.1 x 10 15 cm -2 in that region jRemiian et"al]|2q02b. Finally, 
methy l formate, discovered towards Sgr B2(N) bv iBrown et al] 
has b een also observed in other star formin g regions 
of both high jMacDonald et al] 1 19961 : iGibb et aljboocb and low 
jRemiian & Hollisl l2006t) mass, towards a proto-planetary neb- 
ula jRemiian et al]|2005b. and in comets {Bockelee-Morvan et al] 
l200d : lDespois et al]|200alRemiian et alj|2006b . The column den- 
sity of methyl formate in Sgr B2(N-LMH) is 1.1 x 10 17 cm" 2 
jLiu et alj200lb . 

Sgr B2(N-LMH) is the only source where all three of these 
isomers have been observed. However, while acetic acid and methyl 
formate are concentrate in Sgr B2(N-LMH), glycola ldehyde ap- 
pears to be more diffusely spread through Sgr B2(N) jHollis et al] 
1200 lb. This behaviour is generally shared by other aldehydes 
jSnvdenl2006b . In comets only methyl formate has been observed 
so far. 

An understanding of molecular structure, spectroscopy, and 
photo-absorption processes may be of critical importance in inter- 
preting current observations. In particular, the lifetime of cometary 
molecules versus photo-destruction is a basic parameter for all 
cometary studies: as a matter of fact any error in the photodis- 
sociation rate translates linearly into an error on the abundance 
derived in the cometary nucleus. It is also needed for chemical 
modeling of planetary atmospheres. In this work, we derive elec- 
tronic and vibrational properties of the isomer triplet C2H4O2 de- 
scribed above using the Density Functional Theory (DFT). Photo- 
destruction rates for acetic acid and methyl formate were derived by 
Crovisieri(fl994h using old laboratory absorption data published by 



Suto et all 1 1988b . while the rate for glycolaldehyde is just an esti- 



mate. In section|2]we present a brief outline of the method, together 
with a description of computational settings and results. Section|3] 
contains the application of these results to cometary photochem- 
istry, while discussion and conclusions are in the last section. 



2 THEORY AND RESULTS 



We used DFT (e.g.. Ijones & Gunnarssonl fl989) for the calcula- 
tion of the equilibrium geometry of the electronic ground state 
and of the vibrational spectrum. We then applied the time- 
depen dent extension of the theory (TD-DFT, Mar ques & Gross! 
l2004b to compute the electronic excited states and the resulting 
photo-absorption spectrum for each molecule. 

To obtain the ground-state optimised geometri es we used the 
quant um chemistry program package TURBOMOLE (Ahlric hs et al] 
2007). Technical details about the specific choice of density func- 
tional and atomic basis set can be found in the Appendix. 

After geometry optimisation, we performed the vibrational 
analysis obtaining energies and intensities of the normal modes of 
vibration in the harmonic approximation. Vibrational transitions for 
fundamental configurations of the isomer triplet are given in the 
Appendix. 

Finally, keeping fixed the ground state geometries obtained 
above, we computed the photo-absorption cross-section for each 
molecule. We used two different implementations of TD-DFT in 
the linear response regime, in conjunction with different represen- 
tations of the Kohn-Sham wavefunctions: 

( i) the real-time pro pagation scheme using a grid in real space 
Yabana & Bertschl jl999h, as imple mented in the OCTOPUS com- 
puter program 1 Marques et alfeOoH) : 



(ii) the frequency -space implementation 

(Bauersch mitt & Alrichsl Tl996h based on the linear combina- 
tion of localised orbitals, as given in TURBOMOLE. 

In the real-time propagation scheme (i) the whole 
photo-absorption cross-section of the molecule, up to the 
far-UV, is obtained at once, which is particularly convenient for 
astrophysical applications. Technical details about the TD-DFT 
formalism implemented in the OCTOPUS program are reported in 
the Appendix; the resulting spectra are displayed in Fig.[T] 

In the most widely used frequency-space implementation ( ii) 
the poles of the linear response function correspond to vertical 
excitation energies and the pole strengths to the corresponding 
oscillator strengths. With this method computational costs scale 
steeply with the number of required transitions; electronic ex- 
citations are thus usually limited to the low-energy part of the 
spectrum. Table Q] shows that both combinations BP-TZVP and 
B3LYP-TZVP (see the Appendix for the nomenclature) provide 
a similarly good agreement between the computed transitions and 
the ex perimental results availa ble for acetic acid at room tempera- 
ture jLimao-Vieira et al . 2006), while BP-TZVP results in a much 
closer agreement with the experimental results for glycolaldehyde 
jKarunanandan et al.l2007b . We therefore chose to use the BP func- 
tional for all of our calculations. The resulting BP-TZVP absorp- 
tion spectra for the three molecules, shown in Fig. [2] are obtained 
as a superposition of Gaussian functions with fixed arbitrary widths 
of 0.8 eV. The kind of calculations we performed only yield the 
positions and intensities of vertical, pure electronic transitions, 
and therefore give no information on band widths. However, in 
the av ailable gas-phase sp ectra of acetic ac id jLimao-Vieira et al] 
12006b and glycolaldehyde (Karunana ndan et al]|2007h bands up to 
~ 10 eV show broad profiles with a full width at half maximum of 
about 0.8 eV, produced by the convolution of unresolved vibronic 
structure and the natural width of the transitions. The list of first 50 
excited states and transition intensities is given in the Appendix. 

The two TD-DFT implementations produce compatible re- 
sults in the low-energy region, i. e. up to about 10 eV, while they 
tend to diverge significantly at higher energies. However, TD-DFT, 
as a method, is known to yield dependable results for individual 
transitions only for excitation energies up to the ionisation energy, 
which indeed is close to ~ 10 eV for all three isomers. 

The real-time real-space implementation, on the other hand, 
has been demonstrated to yield good results for the overall den- 
sity of electronic transitions even at high energies, with the caveat 
that single peaks o f fine structure are meaningless (see e. g. 
Mar ques etal] |2003): at energies > 10 eV only the envelope of 
the spectra calculated by octopus is expected to be accurate, while 
the resolved fine structure is largely due to standing waves in the 
finite simulation box, is strongly dependent on the size of the box 
and is only partly quenched by the absorbing boundary conditions 
we adopted. Since we use these spectra at high energies to compute 
absorption rates in a continuous spectrum (see Sect.|3j, the effect of 
spurious fine structure averages out when integrating over ranges of 
many eVs, making these theoretical spectra quite adequate for their 
intended purpose. 



3 PHOTOCHEMISTRY 

Photo-absorption rates of the three isomers are computed by means 
of the relation 



^(s" 1 ) = / S (E)a(E)dE 

J AE 



(1) 
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Figure 1. Comparison between the photo-absorption cross-sections (in Figure 2. Same as Fig. [T]as obtained for the low-energy part of the spectrum 

Megabarns, 1 Mb = 10 — 18 cm 2 ) of acetic acid (red), methyl formate (blue), using the frequency-space TD-DFT implementation of the TURBOMOLE 

and glycolaldehyde (green), as computed up to 100 eV using the real-time package at the BP/TZVP level of theory. 
TD-DFT implementation in the OCTOPUS program. 



Table 1. Comparison between the computed electronic transitions as ob- 
tained using different exchange-correlation functionals and the experi- 
mental data reported by iLimao-VieTr a et al. 1 200(J for acetic acid and by 
iKarunanandan et al] 120071) for glycolaldehyde. 



B-P/TZVP B3LYP/TZVP EXP 





acetic acid 




5.60 


5.84 


6.09 


6.93 


7.49 


7.22 


8.20 


8.30 


8.15 


8.25 


8.50 


8.35 


9.20 


9.25 


8.82 


10.07 


10.54 


10.29 


Avera 


ge relative error 


(%) 


3.4 


3.2 






glycolaldehyde 




4.49 


5.61 


4.51 



Average relative error (%) 
0.4 24.4 



where, for the sake of co mparison, Spj(E) is the solar spec- 
trum at 1 AU provided by Hue bner et M\ dl992h . and a{E) the 
photo-absorption cross-section for a given molecule. The result- 
ing photo-absorption rates are reported in Table [2] where we also 
show results for a radiatio n density expected in a p hotodissocia- 
tion front near an OB star dDraine & BertoldJ[l996h . We also list 
separately the contribution of the Lya line in the solar spectrum 
to the absorption rates, which is of the orde r of ^15% of th e to- 
tal, in agreement with previous estimates of ICrovisIe n dl994l) . As 
evident in Fig. [2] and Table [4] glycolaldehyde presents a relatively 
strong band near 4.5 eV, at an energy ~2 eV lower than the first 
transitions of comparable intensity in the other two members of the 
triplet (namely ~6.9 eV for acetic acid and ~7.5 eV for methyl for- 
mate). Since the solar spectral distribution decreases very steeply in 
this energy range, the estimated glycolaldehyde photo-absorption 
rate in the solar radiation field is by and large dominated by this 



single band, and results to be more than 100 times larger than those 
of acetic acid and methyl formate. Thi s effect is not present for 
photo-absorption rates estimated for the lDraine & Bertoldil d 19961) 
spectrum, since it is much flatter than the solar one in the UV, and 
thus no single band dominates the absorption for any of the three 
species. TD-DFT calculations are known to provide transition en- 
ergies accurate within about 0.3 eV. To test the robustness of our 
results we shifted the energies of all transitions by ±0.3 eV, and 
checked how this affects the resulting photo-absorption rates. Also, 
since TD-DFT provides no information on the intrinsic width of 
the calculated bands, we assumed for all bands a full-width at half 
maximum of ~0.8 eV, consistent with published expe rimental data 
dLimao-Vieira et al 2006; iKarunanandan et a"i]|2007h . The overall 
variations in photo-absorption rates are within a factor of 2, mean- 
ing that our conclusions are rather firm. This accuracy may appear 
surprising for the estimated absorption rates in the Lya line, which 
would naively be expected to vary very strongly with the calculated 
positions of molecular bands, producing either a very large absorp- 
tion rate if it happens to be close to a strong one or a very small 
one if it falls in a gap in the absorption spectrum; however, the as- 
sumption (based on experimental data) of a 0.8 eV FWHM for all 
bands vastly reduces the dependence of absorption rates on small 
(i. e. ~0.3 eV) variations in band positions. 

In a recen t analy sis of ice composition in comet Hale-Bopp, 
ICrovisieret al. (2004) quote a photodestruction rate (at 1 AU) for 
acetic ac id /3cH-jCOOH = 5.1 x 10~ 5 s _1 . This value was pro- 
vided bv lCrovisietl d 19941) . that also reported a rate /3ch 3 ohco = 
4.7 x 10~ 5 s _1 for methyl formate. The photodestruction rate 
of glycola ldehyde is unknown and assumed to be 1 x 10 -4 s _1 
at 1 AU dCrovisieret at] 12004). Assuming a unit photodestruc- 
tion yield and a factor of 2 indeterminatio n in the c alculations, 
our photodestruction rates are consistent with Crovisier ( 1994) val- 
ues, th at were based on old laboratory data provided bv lSuto et al.l 
(1988). The case of glycolaldehyde is different, since most of 
photo-absorption is produced by the band at 4.5 eV, that appears 
to be too low in energy to provide a unit photodestruction yield. 
As a consequence, photodestruction channels are activated just via 
absorption in the high-energy bands, leading to a photo-absorption 
rate /3ch 2 ohcho ~ 3 x 10~ 5 s" 1 (cf. T abled, much lower than 
the rate assumed in lCrovisier et alj d2004l) . 

Photons absorbed in the lower energy bands of CH2OHCHO 
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Table 2. Photo-absorpti on rates (s —1 ) for the isomer triplet C2H4O 2 in two radiation fields, assuming a band FWHM of ~0.8 eV based on available ex- 
perimental data ( Limao-Vieira et al. 2006; Karunanandan et al. 2007). The first column lis ts the computed absorption rates for the whole solar spectrum at 
1 AU fau ebner et alj|l992h . the second column in the solar Lyct line at 1 AU iHuebner et al. 1992), the third in a photodissociation front near an OB star 
iDraine & Bertoldilll996l) . 



Species Solar flux at 1 AU Solar Lya at 1 AU Unshielded OB field 

(Huebner et al. 1992) (Draine & Bertoldi 1996) 

Acetic acid CH3COOH 6.7(-5)W 9.7(-6) 2.9(-9) 

Glycolaldehyde< 2 ) HCOCH 2 OH 9.6(-3) — 1.3(-10) 

2.8(-5) 3.7(-6) 2.9(-9) 

Methyl formate HCOOCH3 5.0(-5) 7.4(-6) 3.0(-9) 

f 1 ' 1.8(-5) = 1.8 x 10~ 6 . 

( 2 ' First row: photo-absorption rates due to the first two absorption bands (up to ~4.5 eV); second row: photo-absorption rates due to the remaining bands. 
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Figure 3. IR emission spectrum of one glycolaldehyde molecule at 1 AU 
from the quiet sun, estimated by assuming that all the energy absorbed in 
electronic transitions is converted in vibrational excitation and subsequently 
reradiated (see text). Emission in each vibrational band is represented by a 
vertical bar whose abscissa corresponds to the wavelength of the vibrational 
mode and whose height equals the calculated emission intensity. 



may be resonantly scattered, or, if the molecule undergoes inter- 
nal conversion, re-emitted in the IR range. Using an adapted ver- 
sion of a Monte Carlo model developed for emission by polycyclyc 
aromatic hydrocarbons, which is indeed assumed to be pumped 
by complete internal conversion of th e energy abso rbed via elec- 
tronic transitions in the visible and UV dMula s 1998), we have con- 
structed the expected IR emission by glycolaldehyde, powered by 
the solar flux at 1 AU, assuming all the energy absorbed in the band 
at ~4.5 eV to be emitted in the IR. The emission coefficient (per 
molecule) is reported in Fig. [5] 

We have thus far considered photo-absorption rates in the 
standard solar radiation field at 1 AU (for the "quiet" Sun), which 
are relevant for the photochemistry of glycolaldehyde, acetic acid 
and methyl formate in our present solar system. However, from the 
study of stellar proxies for the Sun it appears that young solar type 
stars emit high energy photons at a level three to four orders of mag- 



nitude higher than the present-day Sun, both during the pre-main 
sequence phase when the emission is dominated by intense daily 
or weekly flares (Favata et al. 2005), and during the first phases of 
the main sequence (Micela 2002). Therefore, chemical evolution 
can only be understood within the context of the evolving stellar 
radiation environment. 

With out addressing th e problem of molecular survival in a 
disk (e. g. lVisser et alj |2007). we estimate the effect of the extreme 
UV emission (roughly the spectral range between 13 and 100 eV) 
from solar-type stars of different ages, ex ploiting the emission o f 
six stars from the Sun in Time program iGuinan & Ribasll2002l) . 
whose flu xes are assumed to desc ribe the evo l ution of the Sun's 
emission (Ribas et al ] |2005l) . Since lRibasetai]d2005h reported in- 
tegrated fluxes (at 1 AU), the rate coefficents are approximated by 

/3(t)(s- 1 )=E<-)-^|^ (2) 

AE 

where (<t)ab is the average over the energy range AE of the 
photo-abso rption cross-sectio ns shown in Fig.[T] AE intervals are 
taken from Riba s et al] d2005h . Table 4. Beyond 100 eV a collec- 
tive description of the molecule is not anymore necessary, since 
X-ray absorption cross-sections can be closely approximated by 
adding the atomi c cross-sections of individu al atoms bound in the 
molecule (e.g.. lCecchi-Pestellini et af] 2006). We therefore do not 
perform calculations in the X-ray energy range here. Results are 
reported in Fig. [4] It is evident that the high energy tail of stel- 
lar spectrum provides a significant enhancing of photo-absorption 
rates for acetic acid and methyl formate. Glycolaldehyde photon 
absorption rates are not changing too much, although the extreme 
UV contribution to photo-absorption is comparable to absorption 
in the near UV band at ~4.5 eV. In general, due to the increase in 
the stellar high energy component, the photo-destruction rates of 
the isomer triplet members increase roughly two order of magni- 
tude in the environment of a young solar-like star. 



4 DISCUSSION AND CONCLUSIONS 

In this work we studied the photo-physics of the acetic acid 
and its isomers, glycolaldehyde and methyl formate. Computed 
photo-absorption rates are consistent with literature data for acetic 
acid and methyl formate. In the case of glycolaldehyde, for which 
the photo-absorption was completely missing up to the far-UV, 
our calculations indicate that photodestruction is slower than in 
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Table 3. Total energies for the three species considered. The second column 
shows the energy difference S in comparison with glycolaldehyde. 



Figure 4. Photo-absorption rates in the environments of solar-type stars of 
different ages. Top to bottom: glycolaldehyde, acetic acid, methyl formate, 
and glycolaldehyde with the band at 4.5 eV supressed. 



the other two members of the isomer triplet. However, the overall 
photo-absorption rate is much larger, and likely to produce either 
resonant scattering or IR emission powered by near UV solar pho- 
tons. 

Sgr B2 observations of glycolaldehyde show that, unlike 
acetic acid and methyl formate, its emission is extended in the 
surrounding molecular cloud. Such a behaviour, that has not been 
understood so far (e.g. Chengalur & Kanekar 2003), is part of a 
more general problem involving differentiation in isomers, such as 
e.g., isocyanide isomers (CH3CN, CH3NC). The be haviour of this 
isome r triplet has also recently been discussed by Lattel ais"et al.l 
(2009) as a notable exception to what was called the "Mini- 
mum Energy Principle" (MEP) whereby whenever several iso- 
mers are possible for a given formula, their observed abundances 
are in order of binding energy. The authors made the educated 
guess that this anomaly is first created by differences in the 
chemical pathways leading to the formation of glycolhaldehyde, 
acetic acid and methyl formate on the surfaces of dust grains, 
and then preserved due to large energy barriers for the con- 
version among them. Our results may help to shed some more 
light on this problem. Although both hot core and the embedding 
molecular cloud are dark regions, cosmic-rays provide a source 
of UV photons at high visual extinctions by exciting molecu- 
lar hy drogen in the Lyman and Werner bands 'Prasad & Tarafdan 



hy c 

Il983h . This locally generated photon flux typically has fluences 



lower than 10.000 photons cm 2 s 1 dCecchi-P estellini & Aiello 



1992), and may produce imp ortant chemical effects dGredel et al 
1989; Bennett & Raised 120070 . However, the similarity in both 
the photo-absorption cross-sections and ionisation potentials of 
the three species makes chemical differentiation due to selective 
photo-destruction unlikely. As a possible explanation of the ex- 
tended spatial scale of glycolhaldehyde, we consider the possi- 
bility of a slow, selective isomerisation, i.e. the possibility that a 
species may convert itself into another member of the C2H4O2 
triplet by interacting with the radiation field. This would imply 
an isomerisation mechanism which operates on a timescale which 
is much lo nger than the typi cal lifetime of a hot core (few times 
10 4 years. Iwilner et al.ll200ll) but still short enough to be effective 
on the timescale of th e lifetime of a molecular cloud (~ 10 8 years, 
Willi ams etall fcoOO). In this framework, hot cores would reflect 
the relative abundances among the isomers as created by their pro- 
duction mechanisms, whereas in the molecular cloud the isomers 





Total energy (eV) 


<5(eV) 


Glycolaldehyde 


-6231.83 


/ 


Metyl formate 


-6232.32 


0.49 


Acetic acid 


-6233.04 


1.21 



would "relax" to the most stable one, namely glycolaldehyde, in 
agreement with the MEP. 

Isomerisation may be induced by the absorption of radia- 
tion essentially in two main ways; following the absorption of a 
UV-visible photon the molecule could 

( i) move to an electronic state whose energy surface presents a 
minimum close to the equilibrium configuration of another isomer; 

( iij convert, via one or more non-radiative transitions, a substan- 
tial part of the electronic excitation energy into vibrational energy 
allowing the overcoming of the isomerisation potential barriers. 

The present data do not allow a discrimination between the two 
cases. Moreover, the analysis of the first isomerisation chan- 
nel would require a detailed study of the energy hypersurfaces 
in the excited sta t es acc essible with photons generated by the 
IPrasad & Tarafdar] dl983l) mechanism. Therefore, we will discuss 
qualitatively the second isomerisation process. We assume that ev- 
ery absorption heats up the molecule in a time scale character- 
istic of electronic transitions (~ 10~ 8 s), that then decades via 
non-radiative transitions (< 10 -10 s). The electronic excitation 
energy can be uniquely released by a cascade of vibrational transi- 
tions. We also assume, as simplifying hypotheses, that (i) all the ex- 
citation energy is converted in vibrational excitation, (ii) that there 
exist only one isomerisation barrier, (Hi) that the excitation energy 
is far greater than this barrier and (iv) that the isomerisation rate 
is high above the threshold and zero below it. Then every time a 
molecule absorbs an UV-visible photon, it will be vibrationally 
heated. If the energy of the absorbed photon is above the isomerisa- 
tion threshold, the molecule will establish a statistical equilibrium, 
in which the probability to find it in one of the isomeric configura- 
tions will be proportional to the density of vibrational states at the 
given energy. The molecule will then cool down, in timescales of 
the order of a second, with a vibrational cascade (all vibrational 
modes are IR-active for all three isomers). Its energy will thus 
eventually fall below the barrier for the isomerisation. When this 
happens, the proportion among the isomers is frozen because the 
conversion rate among the species drops to zero, the abundances 
are therefore those given by the ratios among the densities of vibra- 
tional states. Since we are here dealing with three conformations, 
one should consider at least three different isomerisation channels, 
each with its different barrier(s). Whatever they are, however, as 
long as they are easily overcome with the energy of a single UV 
photon, the ratio of the abundances of the isomers, if they are ex- 
posed to UV light, should be fixed at the ratios of the densities of 
states at the threshold(s). 

The densities of vibrational states for the three species, 
calculated in harmo n ic app roximation through the algorithm of 
IStein & Rabinovitch] d 1973b . is shown in Fig. \5\ The distances 
among the curves are due to the differences in the total energies 
shown in Table [3] and the differences in the frequencies of the vi- 
brational modes. 
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Figure 5. Densities of the vibrational states for the three species as a func- 
tion of energy. 



The naive model drawn above would essentially turn all iso- 
mers into glycolaldehyde as soon as they absorb an UV photon. 
However, UV absorption, in the weak radiation field produced by 
to cosmic-ray induced H2 fluorescence, occurs on timescales of the 
order of > 50.000 years, that are comparable to the lifetime of hot 
cores. As consequence, MEP should not operate effectively in the 
hot core phase. 

Therefore, one of the simplistic assumption in the naive model 
must be incorrect, slowing down isomerisation. Either internal con- 
version has a low quantum yield for these molecules, meaning that 
they get vibrationally heated only a small fraction of the times they 
absorb an UV photon; or isomerisation is not fast whenever vibra- 
tional excitation is sufficient to overcome reaction barriers, due to 
the morphology of the molecular potential energy surface. The lat- 
ter can easily be the case if the unimolecular reaction path(s) to 
isomerisation are narrow and complicated in terms of phase space 
of the ions. 

Of course, conversion among isomers needs not to proceed 
only via unimolecular reactions induced by radiation: any conver- 
sion reaction including chemical reactions with abundant enough 
partners and without activation barriers could equally well oper- 
ate on the right timescales to fulfil the MEP in molecular clouds 
but not in hot cores. A thorough study of the potential energy sur- 
face of these three isomers, including the reaction paths connecting 
them, is called for in order to further progress in the understanding 
of this observational puzzle. 
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We performed all calculations in the framework of Density Func- 
tional Theory, using t he quantum che mistry packages TUR BOMOLE 
jAhlrichs etal]|2007l) and OCTOPUS jMarques et al]|2003l) . 



4.1 Geometry optimizations 

To find the electronic ground-state geometry of the three molecules 
considered, we first tested different combinations of density func- 
tional and basis set exploiting TURBOMOLE, to identify the most 
suitable one for our purposes. To this aim, we c ompared the experi- 
menta l ground state geometry of the acetic acid dLimao-Vieir a et all 
2006) with the one we obtained using the combinations: BP-SV(P), 
BP-TZVP, B3LYP-SV(P), and B3LYP-TZVP (see the TUR- 
BOMOLE manuaf] and references therein for precise defini- 
tions of basis sets and functionals). All combinations are in 
good agreement wi t h exp erimental and previous theoretical data 
dLimao-Vieira et al.l |2006J), with the hybrid functional B3LYP 
showing a very slightly lower relative average error. This is not 
unexpected, since this functional is known to produce good result s 
with other classes of organic molecules (e.g. jMartin et al.lll996t) . 
Considering that both functionals show the same relative average 
error with both basis sets, we chose the larger TZVP basis set, since 
this is suggested by the TURBOMOLE manual as the default to get 
reliable quantitative results. 



4.2 Vibrational properties 

For our modelling purposes, we calculated the harmonic vibrational 
frequencies of acetic acid, glycolaldehyde, and methyl formate. 
All calculations were performed at the BP/TZVP level using the 
TURBOMOLE program pa ckage, and resulted compatible with pre- 
viously published results dLimao-Vieira et al, 2006; S enen 3 12004 
ISenentetal.ll2005h . The absolute intensities S of the IR-active 
modes are given in units of km-mol -1 . For each vibrational mode 
of frequency v expressed in cm -1 the corresponding Einstein A 
coefficients for spontaneous emission can be computed as: 



A(s~ 



1.2512 x 10" 



8tt .2 
Nag- 



s' 



mol • cm 



cm '/ \ km -mol 1 J km ■ s 
Na being the Avogadro's constant and c the velocity of light. 

4.3 Electronic spectra 

In the real-time implementation of TD-DFT as given in OCTOPUS, 
the time-dependent Kohn-Sham equations are directly integrated 
in real time and the wavefunctions are represented by their discre- 
tised values on a spatial grid. The static Kohn-Sham wavefunctions 
are perturbed by an impulsive electric field and propagated for a 
given finite time interval. In this way, all of the frequencies of the 
system are excited. The whole absolute absorption cross-section 
o-(E) then follows from the dynamical polarisability a(E), which 
is related to the Fourier transform of the time-dependent dipole 
moment of the molecule. The relation is: 

o{E) = ^ (3) 
1 www.turbomole.com 
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where h is Planck's constant, 9{q(_E)} is the imaginary part of the 
dynamical polarisability, and c the velocity of light i n vacuum. 

We perfo rmed the OCTO PUS calculations using iBeckd Jl988h 
exchange and iPerdewl d 198q) correlation functionals. The ionic 
potentials are replaced b y norm-conserving pseudo-potentials 
dTroul lier & Ma rtinslll99lh . We used a grid spacing of 0.12 A and 
a box size of 6.8 A, determined by convergence tests on ground 
state properties and on the photo-absorption spectrum at energies 
< 10 eV. This box size ensures that each atom is at least 4 A 
away from its edges. We furthermore added a 1 A thick absorb- 
ing boundary, which partially quenches spurious resonances due 
to standing waves in the finite simulation box used to c onfine the 
molecules jYabana & Bertscrll999hrMaraues et al.l2003h . We used 
a total time integration length T=20 ft/eV, corresponding to an en- 
ergy resolution of ft/T=0.05 eV. For the numerical integration of 
the time evolution we used a time step of 0.0008 ft/eV, which en- 
sured energy conservation with good accuracy, within numerical 
noise. 

In the most widely used frequency-space TD-DFT implemen- 
tation, the poles of the linear response function correspond to ver- 
tical excitation energies and the pole strengths to the correspond- 
ing oscillator strengths. With this method computational costs scale 
steeply with the number of required transitions and electronic ex- 
citations are thus usually limited to the low-energy part of the 
spectrum. The frequency-space TD-DFT calculations with TUR- 
BOMOLE were performed at the BP/TZVP level of theory, since 
this showed the best agreement with experimental data. We report 
in following Table [4] the first 60 singlet-singlet electronic transi- 
tions of the three molecules under study with the corresponding 
oscillator strengths. 
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Table 4. Vertical electronic transitions (eV) and corresponding oscillator 
strengths of the three molecules considered, as computed at the BP/TZVP 
level using the TURBOMOLE program package. 
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